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Abstract 

In this paper we present a phenomenological treatment of charm quark energy loss before 
fragmenting into D mesons and calculate nuclear modification factor, ’ Raa and azimuthal 
anisotropy, ’tV of D mesons for lead on lead collision at LHC energy of ^fs=2.IQ A TeV. 


1 Introduction 

Whenever two heavy ions collide with ultra-relativistic energies and pass through each other, there 
remains a region at their point of collision where we have a system of partons, composed of gluons, 
and some quarks and anti-quarks in apparently deconfined state but within a very small volume 
in space which is larger than the average nucleon sizes. These quarks and gluons soon relegate to 
their lower energy states via elastic and inelastic collisions and reach an equilibrium momentum 
distribution. This locally thermal equilibrated system is known as quark-gluon-plasma (qgp) [T|— 
[3]. Investigating the properties of quark gluon plasma remains a major activity of present day 
high energy nuclear physics and promises a deeper understanding of the laws of quantum chromo¬ 
dynamics (QCD). Suppression of hadrons [J|, heavy quarks [5], jet-quenching [6], radiation of 
thermal photons [7|, suppression and regeneration of «//?/> {8j are some of the observables which we 
can evaluate as the signatures for the formation of qgp. To start with our discussion let us look 
briefly into the heavy quark production in relativistic heavy ion collision. 

The production rate for the massless gluons and lighter quarks can be traced throughout heavy 
ion collision phases as a minimal momentum transfer is required for their production. On the other 
hand owing to large mass, charm quarks are believed to be mostly produced in the pre-equilibrium 
phase of relativistic heavy ion collision , where partonic momenta are relatively very high. The rate 
of production of charm quark is limited in the later phases of collision history [9], as temperature 
of the thermal medium is far below the charm quark mass. Thus being separated from the bulk 
of qgp, and due to its small numbers, heavy quarks can serve as probes to qgp properties. As 
expected, heavy quarks like other probe particles undergo collisional and radiative energy loss due 
to scattering with the medium partons and this is observed as a suppression in final heavy meson 
spectrum via nuclear modification factor, 1 Raa ’ mm- Evidently recent data from RHIC and 
LHC |13l [Ti j experiments show large suppression for heavy mesons which is almost identical to 
lighter mesons and against the general belief that heavy quark may not loose energy comparable 
to gluons and light quarks. Thus studying heavy quark energy loss mechanisms in thermal bath 
and calculation of nuclear modification factor, , Raa\ have emerged as contemporary topics and 
will also be featured in the present literature. 

Next let us look into another observable, heavy quark azimuthal anisotropy, V 2 also featured in 
this paper. Data observed for non-central collisions both at RHIC and LHC, exhibits considerable 
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elliptic flow particularly for the low momentum D mesons m , and is found to be comparable to 
lighter mesons, V 2 - However one must be careful not to mix heavy mesons results with medium 
partons’ elliptic flow unless heavy quark is assumed thermalized and cannot be distinguished from 
bulk medium any longer. This is however contrary to the conclusions that heavy quark may not 
fully thermalize in quark gluon plasma unless it interacts with the medium strongly m- On the 
other hand with the medium temperatures reaching 400-500 MeV at LHC collider energies, it is 
predicted that low momentum charms may thermalize within the given life time of qgp, which may 
explain D mesons’ large elliptic flow at low transverse momentum sector. 

Thus along with nuclear suppression and elliptic flow, the study of heavy quark evolution in 
quark gluon plasma presents an unique and contemporary field of research. There have been some 
recent theoretical and phenomenological studies of heavy quark dynamics in the thermal medium 
most notably (Em- 121 ]), where along with Raa, and V 2 , other observables such as heavy quark 
correlation, C(A(j),pT,y,r l 7 ) in transverse momentum, azimuthal angle and rapiditites [22], heavy 
quark energy loss per unit path-length in medium, dE/dx, and transport coefficients, such as 
momentum broadening, q, drag coefficient, Ai(p) and diffusion coefficient, Bij(p), heavy quark 
thermalization rates have been calculated in detail. The temperature, momentum and the path 
length dependence of the above transport coefficients using transport models, at both leading order 
and next-to leading order feynman diagrams using resummed loop diagrams calculations in hard 
thermal loop approximation has been calculated, to get a more rigorous picture of heavy quark 
evolution dynamics in quark gluon plasma. 

In this paper we present a phenomenological treatment of charm quark energy loss and calculate 
nuclear modification factor, ’ Raa ’ and azimuthal anisotropy, ’rV of D mesons, at LHC y/s= 2.76 A 
TeV. The recent data from ALICE has shown D mesons, T 2 ’ and ’ Raa ’ f° r various centrality classes. 
This paper tends to explain the data using collisional and radiative energy loss models. Over the 
few following sections, we would discuss about charm production, charm energy loss, calculations 
of the above mentioned observables, followed by discussions of our results and conclusions. 


2 Charm Production 


As discussed in the previous section, charm quark is produced mostly in the pre-equilibrium phase 
of heavy ion collision and its production cross-section can be approximately scaled from proton 
on proton collision using nuclear overlap function, ’ Taa^Y , where commonly called impact 
parameter, is the separation of the centres of two colliding nuclei and denotes centrality class of 
the collision. One can now calculate charm production at the leading order(LO), in , AA > collision 
as 0, 2U, 
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where px and yi .2 are the momenta and rapidities of produced charm and anti-charm and x a and 
x b are the fractions of the momenta carried by the interacting partons from their respective parent 
hadron/nucleon. fi/j(x,Q 2 ) are the partonic distribution functions inside the nucleons. We have 
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Figure 1: Nuclear modification factor, ’i2AA’)(top left)0-10% centr., (top right)0-20% centr. and 
(bottom)30-50% centr. 
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used CTEQ5M partonic structure function [25] with EKS98 shadowing function [26] for our parton 
distribution function, and a constant strong coupling, a s ~0.3 for our calculations. 

We have used TAA{b)= 277 fm~ 2 for 0-10% centrality, 260 fm~ 2 for 0-20% centrality, and 40 
frn -2 for 30-50% centrality for Pb+Pb at LHC, y/s = 2.76 TeV/nucleon calculated from Glauber 
formalism. 

The pt distribution of charm production is calculated using Eq.[2]at the leading order(LO) and 
a Ji-factor of 2.5 is taken to include the next-to-leading order effects. The fundamental processes 
included for LO calculations of invariant scattering amplitudes [29], |M| 2 are g + g —>• c + c;q + q —>• 
c + c and finally the differential scattering cross-section to be used in Eq. [T] is given by, 


dd 

dt 


1 

167t(s — ?n 2 ) 2 


|M| 2 


(3) 


We can now move over to calculation of the energy loss mechanism of the produced charm 
quarks. 


3 Medium effects on charm 

3.1 Collisional energy loss 

The charm quark on entering quark gluon plasma undergoes collision with the medium partons 
with the rate, F ~ j ~ & 2 ^ 2 Pqgp- ’A’ is the mean free path between two successive collision, 
’<72-> 2 ’ is the total elastic scattering cross-section and ' p qg p’ is the qgp medium density, respectively. 
The scattering processes in our calculations involve a small momentum transfer, ’q^’ that slows 
down the charm quark. Thus eikonal approximation for the probe charm quark energy, E q± 
is maintained throughout and 2 —2 scattering processes are calculated within the framework of 
pQCD assuming a constant value of 0.3 for strong coupling, , a s \ 

The assumptions taken above may not be true for all charm energies particularly for low energy 
sector where non-perturbative effects may play a vital role. Similarly higher order corrections and 
loop diagrams resummation techniques also play important roles as demonstrated earlier by m ■ 
We will come back to this later. 

One can approximtely show collisional energy loss per unit length traveled by charm quark in 
the medium as, 


dE co u 
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where n* are thermal medium parton distribution with 3-momentum, E and E' are the energies 
of charm before and after scattering, and ’e = E — E h is the energy loss in a collision. 

The above expression can be simplified using techniques (see ref. [28]) assuming that E, E' T, 
where T is the temperature of the medium. The final integration looks like; 


dE co u 

dx 


-Ww. 
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(5) 


One can now calculate analytically to get the expression by A. Peshier et al., [28] to obtain 
dEcii/ dx. 

The leading-order(LO) scattering processes considered here for differential scattering cross- 
section, dd/dt, are gc —> gc and q(q)c —>• q{q)c and can be calculated or can be taken from 


4 







3.2 Radiative energy loss 

As discussed in the previous section, the charm quark upon entering the hot and dense medium 
undergoes a series of scattering before coming out of the medium as D meson. Upon each of the 
scattering the charm may acquire some virtuality which it looses via radiation of gluons. Radiative 
loss has emerged as one of the prominent energy loss mechanisms of probe charm quark. Radiative 
energy loss along with collisional loss is able to explain experimental data satisfactorily. However 
one must be careful while comparing with various models available for heavy quark radiation 
mechanisms as each of the models has its own set of cuts and constraints imposed in their respective 
calculations (see [301I3T] ). 

Now let us return to the present calculation. Each radiated gluon has its formation time |32] 
after which it is separated from the radiating probe particle and is given by 

2a; 

Tf= 


If V/’ becomes larger than the , ’A’, which is the mean free path or the average distance between 
two successive collisions, then the radiations from successive scattering interfere destructively com¬ 
monly called Landau-Pomeranchuk-Migdal effect |33j . As the result the gluon spectrum from single 
scattering gets suppressed. Following the earlier developments by (see Ref. [3|]) along with the 

hierarchy E,E> 3> q±_ > u 2> k± > /id, I have heuristically assumed the gluon spectrum from 

_ T f _ 


single scattering to be suppressed by no. of coherent centres, , N co f l = . 
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where 
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is the simple Gunion-Berstch formula [35] for gluon radiation from single scattering centre, and 
’R(r/)’ is the dead-cone factor taken from [36], u and r? are the energy and rapidity of the emitted 
gluon, ’L’ is the average path length traveled by charm and , q 1 is the average momentum broadening 
which can be calculated as 



(9) 


and is roughly found to be of the ’0(lGeV 2 /fm)’ [201IM112T1 [23] . C A is the colour factor taken to be 
N c =3.0. The coupling strength, ’a s ’ is taken a constant value of 0.3 throughout our calculations. 
Now the radiative energy loss per unit length traveled by the charm can be written naively, 


dE ra< i (a;) ^ 


dx (A) ’ 

Where the average radiated energy by charm is, 
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And the average mean free path can be written as 
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where 
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Putting the above equations together and after a brief simplification one may finally write eqn. 
(9) as , 
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The limits of the integration for the scattering processes above are calculated to show; 
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Next we can assume as first approximation that the ’total energy loss per unit path-length’ is 
simple addition of collisional and radiative terms shown above and is given by, 


dEi ot dE co n dE ra d 


dx 


dx 


dx 


(16) 


The average total energy loss, ’Aand total momentum loss, ’A pr\ are then calculated over 
average path length, (L), traveled by charm in the transverse direction [24] and can be shown as; 


= f drr J d<j>L{r, 4 >)T aa {t, b ) 

/ drr f dcj>T AA (r, b) 

L(r, cj)) = —r cos cf> + \JR 2 — r 2 sin 2 cj) (17) 

where r , cj) are the position and angle of motion for the charm inside the medium on the transverse 
plane. 

We have assumed only temperature variation of the qgp like medium assuming Bjorken longi¬ 
tudinal expansion only [57] . 

TqTq = T 3 t = TfT f (18) 

No transverse expansion of the medium considered in the present calculation which we feel may 
not be of much effect when a very high momentum probe particle is considered and the calculation 
is done for charm rapidity, |y| ~ 0.0 considering the medium fluid to be always on local rest 
frame [38, 39]. 

The density of the system is roughly assumed to be, 

1 dJV 

"» (T) = <19> 

The charged particle multiplicity, 'dN/dy is taken from [0Oj and multiplied by factors 1.5 get 
dNg/dy. 


6 

























Figure 2: Azimuthal anisotropy, , V 2 i , (left)0-10% centr. and (right)30-50% centr. 


After the energy/momentum loss, we can fragment the charm momentum both from A A and pp 
collisions into D-mesons as D mesons data are readily verifiable from experiments. Schematically, 
this can be shown as, 


d 3 cr d 3 a 

d 3 p Q d 3 pQ 


<8> D(Q 


Hm) , 


( 20 ) 


where the fragmentation of the heavy quark Q into the heavy-meson Hq is described by the function 
D(z). We have assumed that the shape of D(z), where z = Pd/Pc , is identical for all the D-mesons 
and is given by; 



n D 

z[l 1A - W(i - *)] 2 ’ 


( 21 ) 


e p is the Peterson parameter ~0.13 and is taken from m- The normalization condition satisfied 
by the fragmentation function is; 

[ dzD(z) = 1. (22) 

Jo 

With these energy loss mechanisms in hand, we move over to next section calculating charm 
quark observables. 
sectionObservables 


3.3 Nuclear Modification Factor 

The energy loss of quarks (light and heavy) and gluons is easily observed via the suppressed pro¬ 
duction of hadrons measured using nuclear modification factor, Raa- 


Raa{pt, y) 


dN A A/d 2 p T dy 
(■ T A A)do p p/d 2 p T dy 


(23) 


where Naa is the hadron production for the nucleus-nucleus system at a given impact parameter, 
Taa is the corresponding nuclear thickness, and a pp is the cross-section for the production of 
hadrons at the coressponding centre of mass energy/mucleon. 


3.4 Azimuthal Anisotropy 

Non-central collisions of identical nuclei will lead to an oval overlap zone, whose length in and out 
of the reaction plain would be different. Thus, charm quarks traversing the QGP in and out of the 
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plain will cover different path lengths and lose different amount of energy. This would lead to an 
azimuthal dependence in the distribution of resulting charmed mesons, whose azimuthal anisotropy 
could be measured in terms of the V 2 coefficient defined by 


V2{pt) 


f d 0 PT dZd0 cos ( 2 ^ 


f dtp 


dN 


PTdpTdcf) 


where </> is the azimuthal angle in transverse momentum plane. 


(24) 


4 Results and discussion 

To begin with the discussion, let me state that the present calculation assumes a simple phe¬ 
nomenological model for implementing LPM like effects in heavy quark radiation. This is different 
from many earlier calculations such as DGLV, GWP, BDMPS-Z, or ASW [3D1 [3T], @23 — 031 where 
rigorous calculations and implementations of Bethe-Heitler(BH) [45] and LPM effects were done. 
Furthermore there are differences in the various kinematical cuts between these formalisms and the 
present calculation. However in the literature (see Ref. [20]) the authors have stated that there will 
be marginal differences between BH and LPM like effects. All these differences and comparisons 
will feature in our future publications. In this paper I intend to present calculations on nuclear 
modification factor, ’ Raa and azimuthal anisotropy, Vy simultaneously for D mesons compared to 
recent experimental data from LHC. Similar attempts have been done recently in works by mm- 

Let us recall that I have taken LO cross-section for the charm production and multiplied it by 
a K-factor of 2.5 for NLO approximation. This factor gets canceled out when taking the ratios 
such as Raa and ty It is known that dependence on NLO processes changes the magnitude and 
shape of the px distribution of charm cross-section beyond px >20 GeV when compared to LO 
spectrum, and may affect the magnitude of Raa and y at high px region. We will be using NLO 
cross-section for charm production in our future publications. We have also used coupling ’a s ’=0.3, 
thermalization temperature, To=450 MeV, and freeze-out temperature, Tf= 165 MeV respectively 
as the some of the parameters. The transverse expansion of the medium has been neglected as a 
first approximation. 

In Fig. |T] we plot three figures of nuclear modification factor, Raa , of D mesons, each corre¬ 
sponding to different centrality classes. The results are shown for Pb on Pb collision at LHC, y/s= 
2.76 TeV/nucleon. The data are taken from a recent ALICE prelim, data paper [T5] on D mesons. 
The effect on Raa due to collisional loss, radiative loss and total loss are shown separately in each 
of the plots. We notice in all the three plots, there is a particular competition between collisional 
and radiative effects on Raa- At px < 4.0 GeV, collisional loss seems to generate more suppression 
while for px > 4.0 GeV, radiation gives more suppression to charm. But neither collisional loss 
nor radiative loss can explain the data individually and in any case a combination of the two is 
needed for proper explanation of the observed results. In this paper a simple addition of colli¬ 
sional and radiative loss is considered although the respective mechanisms belong to two separate 
models. The author feels that if average energy loss as function of probe charm energies or path 
length dependence of the energy loss is calculated then the difference between collisional and ra¬ 
diative mechanisms can be shown more clearly. These features will be addresses in my upcoming 
publication. 

Let us now returns to the results again. In all three plots we may notice that our calculation of 
Raa involving both collisional and radiative energy loss can explain the data within error bars. We 
may also notice that the heavy quarks are suppressed more in the mid -px region 4.0 GeV < px < 
10.0 GeV than rest of the spectrum. Although there is suppression in the high px region, since 





mass of the charm quark is negligible compared to its momentum, it should behave similar to that 
of light quarks and gluons and there it shows a rising trend but more towards saturation much 
below unity. 

Let us now discuss the plots in detail. The plot (top left) of Fig. [T] is our model calculation for 
0-10% centrality collision (0-7.5% centr. for data). We find that the model calculation shows more 
suppression than the experimental points, for 2.0 GeV< px <4.0 GeV and rises steeply for px < 
2.0 GeV where it is unable to explain the data. But the model is able to explain the experimental 
points above 4.0 GeV charm transverse momentum, within the error bars. Our calculations are 
based on pQCD techniques and applied here for charm quark jets. Hence the present calculation 
may not be able explain results in the low px region where non-perturbative techniques may have 
considerable effects. 

In plot (Top right) of Fig. |T] which is for 0-20% centr. class, the model calculation is able 
to explain results very well for entire px range shown. However we feel that for px > 14.0 GeV, 
since the rise is less, it might underestimate the data. This may be due to the fact that we have 
considered very large emitted gluon energy which might have contributed to the more suppression 
than expected or evident from experiments. Another reason behind this apparent mismatch between 
the present calculation and the data may be due to the absence of NLO processes or higher order 
corrections to the scattering amplitudes which may have some effects on the suppression. 

Also the present calculation assumes eikonal apporximation for the heavy quark propagator 
in the medium. If certain non-eikonal approximation [48] is considered, the effects of soft gluon 
bremsstrahlung at various px regime might change. Similarly transverse expansion of the medium 
might also bring in some considerable effects which is absent in the present calculations. All these 
features are currently under study and will be reported in my future publications. 

Also in plot (bottom) of Fig. [fl we find the rise is even lesser and grossly underestimate the 
data around 12.0 GeV< px <14.0 GeV. The reason here as discussed last may be again due to 
large emitted gluon energy. However in all the three plots, the model calculation explains the data 
to certain extent within error bars. While for low transverse momentum region, it fails as other 
non-perturbative effects etc. must be incorporated. 

Now let us move onto Fig. [2] where we have shown two plots of azimuthal anisotropy V 2 of D 
mesons, for two different centrality classes. Let us recall that our calculation only consists of energy 
loss by charm quarks and no transverse expansion of the medium has been included. In plot (left) 
of Fig. [2] we have V 2 of D mesons for 0-10% centr. Since this is similar to the most central collision 
although not identical to head on collisions, we expect a very small elliptic flow to form and also 
a very small variation with charm momentum (almost flat). This is also evident from our model 
calculation. In plot (right) of the same Fig. [21 we have 30-50% centrality where ellipticity is more 
pronounced and we have more prominent V 2 - Our model calculation roughly explains the data but 
also underestimates it. We feel that absence of any transverse expansion scenario in our model is 
the reason behind the difference between calculated results and the observed data. 

5 Summary and conclusion 

In this paper we have made an attempt in explaining D mesons data at LHC, y/s =2.76 TeV/nucleon. 
We have calculated charm px distribution at leading order multiplied by a k factor. We have in¬ 
cluded heavy quark collisional and radiative energy loss and attempted to calculate collisional and 
radiatve charm energy loss per unit path length. The LPM effect for the radiative processes is 
included heuristically. Consequently Raa and V 2 are calculated for different centrality classes and 
compared with some of the recent D mesons experimental results from ALICE. In future we will 
attempt to incorporate Non-eikonal approximation, soft gluon bremsstrahlung as well as for p % 0, 
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baryonic chemical potential. We would like to include other thermal models or hydrodynamical 
calculations for an expanding medium along transverse direction to give this calculation a more 
realistic picture. Further extension of this work will also include the calculations of other transport 
coefficients such as momentum broadening, q, drag and diffusion coefficients of charm and their 
momentum, temperature, and path length dependencies. 
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